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Abstract 

We examine the effects of a linearly polarized nonresonant radiative field on the 
dynamics of rotationally inelastic Na"^ + N2 collisions at eV collision energies. 
Our treatment is based on the Fraunhofer model of matter wave scattering and its 
recent extension to collisions in electric fields [arXiv:0804.3318vl]. The nonresonant 
radiative field changes the effective shape of the target molecule by aligning it in the 
space-fixed frame. This markedly alters the differential and integral scattering cross 
sections. As the cross sections can be evaluated for a polarization of the radiative 
field coUinear or perpendicular to the relative velocity vector, the model also offers 
predictions about steric asymmetry of the collisions. 
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1 Introduction 

Reactions between ions and simple molecules have been invoked in the chem- 
istry of comets |T], dense interstellar clouds [2], as well as in the atmospheric 
chemistry of planet-like objects, such as lo [3] and Titan Since the late 
1990s, the Na^ -|- X collisions (with X an atmospheric ligand) have been 
recognized to be responsible for the formation of sporadic sodium layers in 
Earth's upper mesosphere [5]. 

The species that take part in the reactions in the upper atmosphere and in 
interstellar space are exposed to electromagnetic radiation of varying inten- 
sity. When polarized, this radiation may create directional molecular states in 
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which the spatial distribution of the molecular axis is itself polarized. This axis 
polarization arises from the nonresonant interaction of the radiation with the 
anisotropic molecular polarizability [6]. Here we examine how such polariza- 
tion may affect the differential and integral cross sections of the ion-molecule 
collisions. 

In particular, we investigate the effect of an intense nonresonant radiative 
field on the rotationally inelastic collisions of Na+ ions with N2 molecules 
at eV collision energies. This collision system is of paramount importance in 
generating sporadic sodium in the mesosphere [5]. 

In our investigation, we make use of a recently developed quantum model of 
collisions in fields [7], which we here adapt for the case of an induced-dipole 
interaction of a nonresonant radiative field with molecular polarizability [B]. 
We limit our considerations to the scattering of molecules by ground-state 
atomic ions. The model is based on Fraunhofer scattering of matter waves [8] , 
[9], [To], and is analytic in both its field-free and field-dependent variant. 

In Section [2| we briefly describe the Fraunhofer model of matter- wave scatter- 
ing and its extension to the case of scattering in nonresonant radiative fields. 
In Section [s} we apply the model to the Na^ + N2 rotationally inelastic 
collisions, and evaluate their differential and integral cross sections and the 
steric asymmetry as a function of the intensity of the radiative field. The main 
conclusions of this work are summarized in Section IH 



2 The Fraunhofer model of matter-wave scattering 

2.1 Field-free scattering 

Inherent to the Fraunhofer model of matter-wave scattering is the energy 
sudden approximation and the assumption of an impenetrable, sharp-edged 
scatterer [7]. As a result, the Fraunhofer amplitude for scattering into an angle 
d from an initial, to a final, |f), state is given by 

/^^fW = (fl/W|i) (1) 

with 

^ j e-'''^^^°"^dR (2) 

the amplitude for Fraunhofer diffraction as observed at a point of radiusvector 
r from the scatterer, see Fig. [Tj Here (p is the polar angle of the radius vector 
R which traces the shape of the scatterer, R = |R|, and /c = |k| with k the 
initial wave vector. Relevant is the shape of the obstacle in the space-fixed 
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XY plane, perpendicular to k, itself directed along the space-fixed Z-axis, cf. 
Fig. HI 



We note that the notion of a sharp-edged scatterer comes close to the rigid- 
shell approximation, widely used in classical [TT], [12], [13], quantum [H], 
and quasi-quantum [15] treatments of field-free molecular collisions, where 
the collision energy by far exceeds the depth of any potential energy well. 

In optics, Fraunhofer (i.e., far- field) diffraction [16] occurs when the Fresnel 
number is small, 

^ = \ < 1 (3) 
rA 

Here a is the dimension of the obstacle, r = |r| is the distance from the 
obstacle to the observer, and A is the wavelength, cf. Fig. [l} Condition ^ 
is well satisfied for nuclear scattering at MeV collision energies as well as for 
molecular collisions at thermal and hyperthermal energies. In the latter case, 
inequality ^ is fulfilled due to the compensation of the larger molecular size a 
by a larger de Broglie wavelength A pertaining to thermal molecular velocities. 

For nearly-circular targets, with a boundary R{(p) = Ro + S{(p) in the XY 
plane, the Fraunhofer integral of Eq. ^ can be evaluated and expanded in a 
power series in the deformation S{ip), 

m = /o(^) + 6) + 6') + --- (4) 

with fo{^) the amplitude for scattering by a disk of radius Rq 

(5) 

and /i the lowest-order anisotropic amplitude. 



2lT Jo 

where Ji is a Bessel function of the first kind. Both Eqs. ([s]) and ^ are appli- 
cable at small values of ^9 < 30°, i.e., within the validity of the approximation 

The scatterer's shape in the space fixed frame, see Fig. [T| is given by 

Ria, P, r, 0,V) = T. ^.u^^MPl)y^piO, ^) (7) 

where (a,/3, 7) are the Euler angles through which the body-fixed frame is 
rotated relative to the space-fixed frame, {6, ip) are the polar and azimuthal 
angles in the space-fixed frame, ^p^,{af3'y) are the Wigner rotation matrices, 
and Ski/ are the Legendre moments describing the scatterer's shape in the 
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body-fixed frame. Clearly, the term with z/ = corresponds to a disk of radius 

i?o - ^ (8) 
V47r 

Since of relevance is the shape of the target in the XY plane, we set ^ = | in 
Eq. ([7]). As a result, 

6{^) = R{a,p,r, f -Ro = R{¥^) -Ro = Es«.^;.M7)>^Kp(f ,<^) (9) 

KUp 

By combining Eqs. ([l]), ([6]), and ^ we finally obtain 

/wf(^) ^ (fl/o + = (f|/i|i) = E H..(f|^;ji)F.,J|,|(A:i?o^) 

■^'^ KUp 

K+p even 

(10) 

where 

' i-mn {-^)h-^ri^l^^ for. + p even and .>p 

F^p= < 

elsewhere 

(11) 

For negative values of p, the factor {—i)'^ is to be replaced by i'^. 



2.2 Scattering in a radiative field 



When subject to an external electric field, the electronic distribution of any 
molecule becomes polarized to some extent. This interaction, governed by the 
molecular polarizability, results in an induced dipole moment. While for the ex- 
perimentally feasible static fields such induced moments are very weak, sizable 
dipole moments can be induced by a radiative field. If the induced-dipole inter- 
action is anisotropic and sufficiently strong, the molecular rotational states un- 
dergo hybridization (coherent linear superposition) which aligns the molecular 
axis along the field vector [6]. The strength of the interaction is characterized 
by a dimensionless parameter Au 

"^^ — -^ <i2> 

with Att = a|| — ax_ the polarizability anisotropy, ol\\^^ the polarizability com- 
ponents parallel and perpendicular to the molecular axis, B the rotational 
constant of the molecule, / the radiation intensity, and e the amplitude of the 
corresponding oscillating electric field. The induced-dipole interaction couples 
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states of the free-rotor basis set with same M but with J's that differ by 0, ±2. 
Thus the resulting hybrid states take the form 



|J,M; Atu) = ^ a;^jv^(Acj)|J,M) for J even (13) 

J=2n 

\J,M;Auj)= «jm(^^)I^.^) for J odd (14) 

J=2ra+1 

where 2n = m+|M| and 2?t,-|-1 = m+|M| with m either 0, 2, 4 . . . or 1, 3, 5 ... . 
The hybridization coefficients aj]^.j{Auj) depend solely on the interaction pa- 
rameter Au. The symbol J denotes the nominal value of J that pertains to 
the field-free rotational state which adiabatically correlates with the hybrid 

state, 

\J,M,Au ^0) ^\J,M) (15) 



Since the hybrid wavefunctions, Eqs. (13) and (14), comprise either even or 
odd J's, the states have definite parity, (—1)"^. 

Apart from possessing a particular energy level pattern, the | J, M, Au) eigen- 
states are aligned along the electric field vector, e. The degree of alignment 
depends on the values of J, M, and Auj. In such states, the molecular axis 
librates about the field direction like a pendulum, and so the hybrid states are 
referred to as pendular. It is the directionality of the pendular states that enters 
the field-dependent Fraunhofer model and distinguishes it from the field-free 
model, which assumes an isotropic distribution of the molecular axes. The di- 
rectional properties of pendular states are exemplified in Fig. |2| which shows 
polar diagrams of both field-free and pendular wave functions at Auj = 25. 

The scattering process in the field consists of the following steps: A molecule 
in a free-rotor state | J, M) enters adiabatically the radiative field where it is 
transformed into a pendular state \J,M,Auj). This pendular state may be 
changed by the collision in the field into another pendular state, \J', M', Auj). 
As the molecule leaves the field, the latter pendular state is adiabatically 
transformed into a free-rotor state \J',M'). Thus the net result is, in general, 
a rotationally inelastic collision, \J,M) — ^ \J',M'). 

In order to be able to apply Eq. ((To|) to collisions in the radiative field, we have 



to transform Eqs. (13) and (14) to the space-fixed frame XYZ. If the electric 
field vector is specified by the Euler angles {ip^, 0^, 0) in the XYZ frame, the 
initial and final pendular states take the form 

|i) = \J,M;Au) = Y.aij,,{Auj)J2^^\,{^,,9,,0)Yj^ie,^) (16) 
J i 

(fl ^ {J',M'-Au\ =J24lA^^)T.^iM'iVe,0s,0)Yj\,ie,^) (17) 

J' 5' 

which is seen to depend solely on the angles 9 and ip. 
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On substituting from Eqs. (16) and (17) into Eq. (10) and its integration, 
we obtain a general expression for the Fraunhofer scattering amplitude in the 
field, 

K+p even 



X ^ aj^j{Auj)bj',lj,{Auj)^m^C{JKj'; 000)C{JkJ'; MAMM') (18) 



jj' 



where AM = M' - M and C( Ji, J2, ^3; Mi, M2, M3) are Clebsch-Gordan co- 
efiients [20]. Since the ion-linear molecule potential is axially symmetric, only 
the coefficients contribute to the scattering amplitude. 



Eq. (|18j) simplifies for special cases. If we limit our considerations to homonu- 
clear diatomics, only the E^q coefficients for even k contribute to the expansion, 
Eq. ([7]), and, consequently, to the scattering amplitude, Eq. ( jlsl . Furthermore, 
if we fix the initial molecular state to the ground state, | J, M) = |0, 0), and re- 
strict the polarization of the radiation in the space-fixed frame to a particular 
geometry, the problem simplifies as follows: 

(i) For a polarization vector collinear with the initial wave vector, s \\ k, 
we have 6'g — > 0, — > 0. As as result, only the p = —AM' term yields a 
nonvanishing contribution and so 



K even 



X E«%(^)^j;M'Hy^^^(^«^'^^000)C(J/cJ^OM'M') (19) 

We see that the angular dependence of the scattering amplitude for the parallel 
case is simple, given by a single Bessel function, J\m'\- 

(ii) If the polarization vector is perpendicular to the initial wave vector, e ± k, 
we have 9^ ^ ^, ^Pe ^ 0. Hence 

Ctj'M'(^) = "1^° E dl,M' (2) ^^oF^MkRo^) 

K,p even ^ ^ ^ 

X ga%(a;)6];;,,(^)^|^C(J«:J';OOO)C(J«:J';OM'M0 (20) 

where d'^^ j^t are the real Wigner rotation matrices. Since the summation mixes 
different Bessel functions (for a range of p's), the angular dependence of the 
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scattering amplitude in the perpendicular case is more involved than in the 
parallel case. 



The Clebsch-Gordan coefficient C{JkJ'; 000) in Eqs. (19) and (20) is nonzero 
only if J + J' is even, since the summation includes only even-K terms. More- 
over, given the definite parity of the pendular states, we see that only parity- 
conserving transitions are allowed, namely J = — > J' = 2, 4, 6, . . . for our 
choice of the initial state. 



We can also see that, for either geometry, only the partial cross sections for the 
J = 0, M = —>■ J' , M' collisions with M' even contribute to the scattering. 
This is particularly clear in the s || k case, where the -FkM' coefficients vanish 
for M' odd. In the £ ± k summation over p arises. Since for k even and 

M' odd the real Wigner matrices obey the relation d'^p j^p (^|^ = —dp^j^, (^|^, 
the sum over p is zero and so are the partial cross sections with M' odd. 



3 Rotationally inelastic collisions of Na+ with N2 in a radiative 
field 



Here we apply the model to the Na"*" + N2(J = J') collisions. The 
polarization anisotropy Aa = 0.93 and rotational constant B = 1.9982 
cm~^ make the N2 molecule a suitable candidate for an experiment on laser- 
assisted ion-molecule collisions. 

According to Ref. [I9], the ground-state Na"'"-N2 potential energy surface has 
a global minimum —2712 cm~^ deep. The effect of this attractive well is neg- 
ligible for low-energy collisions; we chose a collision energy of 5 eV, which 
corresponds to a wave number k = 173.8 A~^. The "hard shell" of the poten- 
tial energy surface at this collision energy is shown in Fig. [3j We found it by a 
fit to Eq. ([7]). The S^o coefficients are listed in Table [l] Due to the sym- 
metry of the potential energy surface, only even-K terms arise. According to 
Eq. (|8]), the Soo coefficient determines the hard-sphere radius Rq, responsible 
for elastic scattering. 



3.1 Differential cross sections 



The field-free state-to-state differential cross section, 

^",Lj',A/'W = l/0,0-.J',M'Wr 



(21) 



see Eq. (10), is proportional to Sj,q, which means that the shape of the repul- 
sive potential provides a direct information about the relative probabilities of 
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the field-free transitions and vice versa. Since for the Na^-N2 system the ^2,0 
coefficient dominates the anisotropic part of the potential, see Table [T| the 
corresponding J = — J' = 2 collisions are expected to dominate the inelas- 
tic cross section. Because of the symmetry, there are no parity-breaking 
J = —>■ odd J' collisions in the Na"'"-N2(J = J') system. 

After averaging over M' and invoking the asymptotic properties of the Bessel 
functions ^21j , we obtain for the parity-conserving J = — >■ even J' collisions 

Xfl^,(^)~cos2(^A;i?o^-0 (22) 



The elastic differential scattering cross section, cf. Eq. (5 ), has a sin 
asymptote, and so is seen to be shifted with respect To the differential cross 
sections for even- J' transitions by a quarter of a wavelength. Known as the 
"Blair phase rule," the shift is a conspicuous feature of Figs. |4] and [5j 

The state-to-state differential cross sections for scattering in a radiative field 
parallel (e || k) and perpendicular (e ± k) to the initial wave vector are given 
by 

^oii''^^W = E^i^i',M'W (23) 

M' 

with 



= /oTj/m'(^) (24) 

The differential cross sections for the Na"*" + N2 collisions are presented in 
Figs. |4] and [5] for an interaction parameter Aa; = 10 and 25, corresponding to 
laser intensities of 2.15 x 10^^ W/cm^ and 5.37 x 10^^ W/cm^, respectively. The 
figures show that a radiative field on the order of 10^^ W/cm^ dramatically 
alters the magnitudes of the differential cross sections, but does not produce 
any "phase shift" of the angular oscillations. Such a "phase shift" is absent 
because only even Bessel functions, which have a cos^ {kRo-d — |^ asymptote. 



contribute to the scattering at any field strength, see Eqs. (19) and (20). 



3.2 Integral cross sections 



The angular range, i) < 30°, where the Fraunhofer approximation applies the 
best, comprises the largest impact-parameter collisions that contribute to the 
scattering the most, see Figs. |4] and [5j Therefore, the integral cross section can 
be obtained to a good approximation by integrating the Fraunhofer differential 



cross section, Eq. (23), over the solid angle siwdd'ddLp, 

c2tv 







^^11;,^)=/ rf<^ / I^:i^^j^\^)smM^ (25) 
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The integral cross-sections thus obtained for the field parallel and perpendic- 
ular to the initial wave vector are presented in Fig. [6j One can see that the 
state-to-state cross section for the J = — >■ J' = 2 collisions steadily decreases 
with the interaction parameter Aa;, whereas the other state-to-state cross sec- 
tions show a non-monotonous dependence. These features can be explained by 
the field dependence of the overlap of the hybridization coe 
and 6j'^,(Au;), affecting the scattering amplitude, cf. Eqs. 
Ref. [Zl. 



icients, a jjj{Aujj 
igl) and Em and 



3.3 Steric asymmetry 



We define the steric asymmetry as 

S„, = (26) 

cry + cr± 

where the cross sections cr||^_L correspond, respectively, to s || k and e ± k. 
The dependence of the steric asymmetry on the induced dipole interaction pa- 
rameter Aa; is presented in Fig. [7j One can see that a particularly pronounced 
asymmetry obtains for the J = —>■ J' = 6 channel. This can be traced to 
the field dependence of the corresponding integral cross sections. Fig. |6| In- 
deed, a conspicuous feature seen in Fig. |6]is the significant dependence of the 
J = — J' = 6 channel on the collision geometry. The integral cross section 
for the J = — s> J' = 6 channel is always greater for the £ ± k geometry 
than it is for s || k because of the non-vanishing d'^p^ii^) Wigner matrices. 



cf. Eqs. (19) and (20). The asymmetry for the J = 0— >J' = 2or4 channels 



is less pronounced, as only terms with M' up to 2 or 4 are involved for e ± k. 

We note that within the Fraunhofer model, elastic collisions do not exhibit 
any steric asymmetry. This follows from the isotropy of the elastic scattering 



amplitude, Eq. (20), which depends on the radius Rq only: a sphere looks the 



same from any direction. 



4 Conclusions 



We made use of the Fraunhofer model of matter wave scattering to treat 
rotationally inelastic ion-molecule collisions in nonresonant radiative fields. In 
accordance with the energy sudden approximation, inherent to the Fraunhofer 
model, the interaction must be dominated by repulsion, which is typically 
well satisfied for ion-molecule collisions down to collision energies on the oder 
of 1 eV. The Fraunhofer model is also inherently quantum and, therefore, 
capable of accounting for interference and other non-classical effects. The effect 
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of the radiative field enters the model via the directional properties of the 
molecular states created by the field. Even a small alignment of the molecules 
was shown to cause a large alteration of the differential and integral cross 
sections. The strength of the analytic model lies in its ability to separate 
dynamical and geometrical effects and to qualitatively explain the resulting 
scattering features. These include the angular oscillations in the state-to-state 
differential cross sections or the rotational-state dependent oscillations in the 
integral cross sections as a function of the intensity of the radiative field. 

We hope that the model will inspire new experimental work based on a com- 
bination of ion traps with molecular beams |22]- 
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Table 1 

Hard shell Legendre moments H«;o, Eq. ([7|, for the Na"^ — N2 potential at a collision 
energy of 5 eV. All odd moments are zero. 





H«o (A) 





6.1221 
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0.5301 
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-0.0359 


6 


0.0022 
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0.0002 
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Fig. 1. Schematic of Fraunhofer diffraction by an impenetrable, sliarp-edged obstacle 
as observed at a point of radius vector r(X, Z) from the obstacle. Relevant is the 
shape of the obstacle in the XY plane, perpendicular to the initial wave vector, k, 
itself directed along the Z-axis of the space-fixed system XYZ. The angle (p is the 
polar angle of the radius vector R which traces the shape of the obstacle in the 
X, Y plane and is the scattering angle. See text. 
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Fig. 2. A comparison of the moduli of the free rotor wavefunctions \J,M = 0), panel 
(a), with the moduli of the pendular wavefunctions \J, M = 0; Auj = 25), panel (b). 
The polarization vector e of the radiative field is also shown. 
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270° 

Fig. 3. Equipotential line R{0) for the Na"*"- N2 potential energy surface at a collision 
energy of 5 eV. The Legendre moments, Eq. ([7|, of the potential energy surface are 
listed in Table [H 
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Fig. 4. Differential cross sections for the Na"*" + N2 (J = — > J') collisions in a 
radiative field for Aa; = 10 (red dashed line) and Aw = 25 (blue solid line), parallel 
to the initial wave vector. The field- free cross sections are shown by the green solid 
line. 
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Fig. 5. Differential cross sections for the Na^ + N2 (J = — > J') collisions in 
a radiative field for Au = 10 (red dashed line) and Auj = 25 (blue solid line), 
perpendicular to the initial wave vector. The field-free cross sections are shown by 
the green solid line. 
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Fig. 6. Partial integral cross sections for Na^ + N2 (J = ^ J') collisions in a 
radiative field parallel, panel (a), and perpendicular, panel (b), to the initial wave 
vector. 
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Fig. 7. Steric asymmetry, as defined by Eq. (26), for Na"*" + N2 (J 
collisions. 
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